Singularities of solutions to Schrodinger equation 
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Abstract 

In this paper we study microlocal singularities of solutions to Schrodinger 
equations on scattering manifolds, i.e., noncompact Riemannian mani- 
folds with asymptotically conic ends. We characterize the wave front set 
of the solutions in terms of the initial condition and the classical scatter- 
ing maps under the nontrapping condition. Our result is closely related to 
a recent work by Hassell and Wunsch, though our model is more general 
and the method, which relies heavily on scattering theoretical ideas, is 
simple and quite different. In particular, we use Egorov-type argument in 
the standard pseudodifferential symbol classes, and avoid using Legendre 
distributions. In the proof, we employ a microlocal smoothing property 
in terms of the radially homogenous wave front set, which is more precise 
than the preceding results. 

1 Introduction 

We consider Schrodinger operators on noncompact manifolds, which is called 
scattering manifolds following R. Melrose (|12j). Scattering manifold is a natural 
generalization of asymptotically Euclidean manifolds, though it is denned as a 
manifold with boundary with Riemannian metric of a special form. 

Let M be a manifold, and we suppose M is the interior of a compact manifold 
M with smooth boundary dM, i.e., M = (M)°, where E° denotes the interior 
of E. Let fibea neighborhood of dM in M that is diffeomorphic to (0, 1] x dM, 
where {1} x dM corresponds to dM. By a map: 



we identify 



= (1-t)-\ (0,1]^[1,oo) 



:= Mnft^ (1, oo) x dM. 
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Such a diffcomorphism is called a boundary decomposition. In this paper we 
fix a boundary decomposition. We also use an atlas of M^ such that each 
coordinate neighborhood is diffeomorphic to (1, oo) x U, where U is a coordinate 
neighborhood of dM. We write 

M = M U Moo, 
where M <s M is a relatively compact open submanifold. 

Definition 1.1 Let g° , g be Riemannian metrics on M. g° on M is called 
conic (with respect to a given boundary decomposition), if there are R > and 
a Riemannian metric h on dM such that 

g° = dr 2 + r 2 h 3k (9)d9^d9 k for (r, 9) G [R, oo) x dM. 

A Riemannian metric g is called scattering metric of long-range type, if we can 
write g = g° + m, where g° is a conic metric, and m is of the form 

m = m°(r, 9)dr 2 + rm){r, 9){drde ] + dO 1 dr) + r 2 m 2 k {r, 6)d6 :i d6 k 

such that for some constants \x\ > and any indices k and a, 

\d k r d%m\r,6)\<C Kka r-^- k for (r, 9) E (1, oo) x K, 1 = 0,1,2. 

Here K is any compact subset of a coordinate neighborhood U C dM. (The 
dependence of constants Cxka on K will be suppressed, when there is no confu- 
sion.) In particular, if /i > 1 and \i\ > 1/2, then g is said to be of short-range 
type. 

In this paper we assume g is a scattering metric of short-range type. 
We consider solutions to the Schrodingcr equation: 

d 

i—u(t,x) = (-A g + V)u(t,x), u(0,-)=u o eH, (1.1) 

where 

I n 

A 5 = — d j9 jk V9~d k , (g lk ) = (g 3k )-\ g = det{g jk ), 

V9 j,k=l 

is the Laplacc-Bcltrami operator, and the function space is given by 

H = {L 2 (M;y/gdx), (•,•)«), {u,v) H =l u(x)v(x)y/g(x)dx. 

Jm 

We suppose V is of smooth short-range type in the sense that V £ C°°(M;R) 
and for some ^ > 1 

\d k d%V(r,9)\<C ka r 2 -^- k . 
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We may assume 



< /JL = (Jl® — 1 = 2/xi - 1 = fl-z = fa — 1 < 1. 

Under these settings the Schrodinger operator H = — A g + V is essentially 
self-adjoint, and we denote the unique self-adjoint extension by the same symbol 
H . Then the equation (jl.ip has the unique solution u(t, ■) = e~ uq for any 
initial data uq e H. 

We now fix (xo, £°) € T*M and to > 0. Our main result concerns a necessary 
and sufficient condition for (&o,£ ) S WF(e~ I *° i? tto) in terms of the scattering 
data for the classical trajectory with the initial condition (x^,^ ). We denote 
by (x(t; xo, £(t; xo, £ )) = exptH p (xo,t; ) the solution to the free Hamilton 
equation on (M,g): 

with initial data (z(0),£(0)) = (xo,£ )• We say (xo,£o) is backward nontrapping 
if x(t;xo,^o) escapes from any compact set in M as t — > — oo, and denote the 
set of all backward trapping points by 7^ . In the next section we show that if 

(aJo,C ) G T * M \ 71 , thcn tnc limits 

r_(x ,e°)- lim (r(i;x ,C°)-2t i0 (<; a ;o^ )), M*o,£°) = Mm 0(t;x o ,f), 

c — > — oo t — > — oo 

p-(xo,€°)= lim p(t;a;o,e ), ^-(zo^ ) = Hm u{t;x ,£°) 

t — > — oo r — > — oo 

exist, where (r,8,p,u>) G T*Moo are local coordinates associated with the 
boundary decomposition (r, #) € M^, 

The limits (r_ ,6-, p~ ,lu~) are the classical scattering data at time — oo. We 
note p~ is the negative square root of the total energy; 9- is the asymptotic 
direction; oj~ is the impact parameter (up to a scaling factor); and r_ is the 
time-shift of the trajectory (times velocity). 

We then set up a reference free system in the quantum mechanics. We set 

Mfree = M X dM, 

Hf Iee = (L 2 (M llcc , VhdrdO), (•, -)free), 



(u,u)free= / u(r, 6)v(r, 9)^h(6)drd9, h = det(hjk), 

JM !ree 

where (hjk) is the Riemannian metric on dM, appearing in Definition 1.1. We 
now set our free Hamiltonian as 

d 2 

H = -g-jj) D(Ho) = {u e H{ree\ H U G Hhee} ■ 

Then H is self-adjoint on Hf lcc . We choose j € C°°((l, oo)) such that 

m = i 1; ifr>2 ' 

J( > 10, if r < 3/2, 



3 



and define J : H{ ICC — > H by 

(J )( ' ~ \ 0, if x^ Moo. 1 J 

Using the natural identification: C -Mf ree , we may write 

(r_,6>_,p~,w~) = lim exp(-tiJ p2 ) o exp tH p (x , £°) G T*M free . 

Now we state the main result of this paper. 

Theorem 1.2 Assume g is a scattering metric of short-range type, and let 
H = — A g + V where V is a smooth short-range type potential. Let uq S TL and 
to > 0, and suppose (&o>£ ) G T*M \ 71. TTien 

Remarks 1. The theorem describe the wave front set of e~ lt ° H uo in terms of 
the classical (inverse) wave operator and e~ lt ° H ° J*uq. Since Hq is essentially 
one-dimensional free Schrodinger operator, we can consider the time evolution 
as a linear transform: (x, £) i— > (x + 2i£, £) in the phase space, and the right 
hand side is easy to characterize. (We can also easily write the integral kernel 
of e - itHo explicitly.) 

2. The theorem suggests that e ltH ° J* e~ ltH is a Fourier integral operator corre- 
sponding to the canonical map: (xo,£o) >- » (f-> P -, This is analogous 
to what Hassell and Wunsch proved in a different setting 0]. We believe that 
this statement can be proved using our method, and it will be discussed in 
forthcoming papers. 

3. As a simplest example, we consider the Euclidean space (M, g) = (M. n ,dx 2 ), 
and H = — A on L 2 (R n ). (R n , dx 2 ) is a scattering manifold of short-range type, 
identified with the interior of the half sphere M = §" = {w € R n+1 ; \w\ = 
1, w n+1 > 0} through the stereographic projection: 

SP:R"-»S£, I MSP(i) = ( 7 ii) 1 (x) = (1 + M 2 ) 1/2 . 

(x) (x) 

Our comparison system is 

ffo = -T7T onL 2 (RxS n -\dr-v^-i), 
ar A 

where Vg»-i is the standard density on S™ _1 and Theorem 11.21 gives a charac- 
terization of WF(e lfA uo) in terms of WF(e l * a T Ii ^j(r)ito). 

4. In [T3] and [T5], one of the authors studied the same problem for the asymp- 
totically Euclidean case, and our result may be considered as a generalization 
of them. However, in [13], Ho = —A is used as the free system, and hence the 
result in [13] is not a special case of Theorem 1.2. The recent paper [15] concerns 
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the long-range case, and the comparison evolution is given by a modified free 
motion. 

The analysis of microlocal singularities of solutions to Schrodinger equations 
with variable coefficients was introduced by a seminal paper by Craig, Kappeler 
and Strauss [1] in 1995. They considered Schrodinger equation with asymp- 
totically flat metric of short-range type on the Euclidean space, and showed 
so-called microlocal smoothing property, that is, the microlocal smoothness fol- 
lows from the rapid decay of the initial state in a conic neighborhood of the 
asymptotic velocity (as t — ► — oo). This result was generalized to Schrodinger 
equations on manifolds with scattering metric by Wunsch [21j , and to long-range 
type perturbation (on the Euclidean space) by Nakamura [2] . We remark that 
both results are more precise than [1] in the sense that the conditions on the 
initial states are time-dependent. Wunsch used quadratic scattering wave front 
set and Nakamura used homogeneous wave front set to formulate the assump- 
tion. The definitions of these notions differ considerably, but it was shown by 
Ito [B] that they are essentially equivalent (up to a linear transform). These 
results are generalized to the microlocal analytic singularities by Robbiano and 
Zuily [HI [UlIIS], and Martinez, Nakamura and Sordoni [TO] . 

The microlocal smoothing properties give us sufficient conditions for the mi- 
crolocal smoothness of solutions, but they do not give us necessary conditions, 
which was addressed by Hassell and Wunsch [3] 0] . They considered the prob- 
lem on manifold with scattering metric of short-range type (somewhat stronger 
assumptions), and gave a necessary and sufficient condition of the microlocal 
regularities of solutions in the nontrapping area. The same problem was con- 
sidered in [13\ 1 1 5 j , and essentially equivalent characterization of the microlocal 
singularities of solutions was proved for Schrodinger equations of long-range type 
on the Euclidean spaces. An analytic analogue of these results (short-range type 
on the Euclidean space) was recently proved by Martinez, Nakamura and Sor- 
doni [11] . The method of this paper is closely related to that of [13j [15] , though 
the geometric structure is quite different. In particular, in [T2], the standard 
classical scattering theory is used to construct the scattering correspondence, 
whereas we essentially construct the classical scattering in the polar coordinate, 
which makes these results being not equivalent in the case of Euclidean space. 
We also note that our method is partially inspired by works on the character- 
ization of the singularities for perturbed harmonic oscillators ([22],[8 ,[2 ). We 
refer pQ, [4] and [15] for other references. 

Compared to the results by Hassell and Wunsch [4] , our result is formulated 
using more elementary terms, and the proof seems simpler. Hassell and Wunsch 
used the multiplication operator e lr l 2t to compensate the oscillation of the 
solution at oo, whereas we use the one-dimensional free Schrodinger evolution 
e -itH f or -j-^at purpose. Thus the correspondence of the microlocal singularities 
is given by the classical scattering relation in our paper, whereas the so-called 
sojourn relation is used in [1]. Also, the scattering wave front set (due to 
Melrose [12]) is used to characterize the wave front set in [4], whereas we use 
the standard wave front set. We note that Hassell and Wunsch constructed a 
parametrix of the evolution operator as a Legendre distribution, whereas we 
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have not used any parametrix construction but an Egorov-type theorem. 

The paper is organized as follows: In Section 2 we study the classical scat- 
tering on manifold with scattering metric. Section 3 is devoted to the proof of 
the microfocal smoothing property of our evolution operator. The smoothing 
property proved in this section is more precise than the previous results, and it 
is of interest in itself (Theorem 3.2), though we need only its very weak version 
in the proof of our main result. We prove the main theorem in Section 4. 

Acknowledgements. SN wishes to thank Andrew Hasscll and Jared Wunsch 
for valuable discussions and comments. He is also grateful to Kenji Yajima for 
his encouragement. 

2 Classical trajectories and scattering data 

In this section we consider the asymptotic behavior of the classical flow gener- 
ated by the kinetic energy function p(x,£) on T*M. Let (M,g) be a scattering 
manifold of short-range type, and (xo,£°) € T*M \ 71, then 



exptH p {x ,Z a ) e T*M X 
for large — t > 0. Thus we can write 

where (r(t; x , 6(t; x , p(t; x , x , £ )) satisfies 




dp 





dp 




dp 



(2.1) 



with 



P(r, 0, p, u>) 



= p 2 + ±v 



(0)ujj(Jk + a (r, 6)p 2 + -a{(r, 6)puj j + — a J 2 k (r, Q)ujjU k - 



By the Cramer's formula we can compute the inverse matrix to 




1 + m 1 

rm 1 



o 




so that we obtain 



\d*d%a (r,8)\ <C ka r- x -»- k , 
\d k r d% ai {r,6)\ < C ka r-( 1+ ^/ 2 - k < C ka r-»- 
\d k d%a 2 (r,d)\<C ka r-»- k . 
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Then (|2.1[) is written explicitly: 



r = 2p + 2a p+ -a{ojj, (2.2) 
9 j = ^h jk w k + l -a? lP + ^a{ k u k , (2.3) 



(2.4) 



2 -u dag 2 1 7 

1 da{ 2 jk 1 <9ai fe 

1 <% feZ <9a 2 lda\ 1 <9< 

- W * w i--5S7P ( 2 - 5 ) 



J r 2 00J afl» ' r 863 ' K r 2 dO^ 

Lemma 2.1 If (x , £°) 6 T*M\T_, then there is large C > smc/i i/iai 

r(t;a;o ) f )>C- 1 |t|-C fort<-C. (2.6) 

T/ie estimate holds locally uniformly in the initial data (xo,£,°). In particular, 
T*M \T- is open. 

Proof. First note that 

p{r(t; xo, 0(t; z , p(«; so, w(*5 ^o, £ )) = p(*o, £°) = Po 
by the conservation of the energy. Then, since 

Po = {\ + a )p 2 + ^(h jk - ^a{a k + af)w,-w fc + \{p+ -afci? 

> (g + fl o)-° 2 + ^ - ^«i a i + 4 k )vjVk, 

we obtain 

|p(*;z ,£ )| < C, k(t;^ ,C°)| < CY(i;z ,£°) (2.7) 
as long as r(i;xo,£°) is large. By direct computation we have 



;= 21, f 2^rp + 2^¥ p + 2a p - ^a{c 



1 . 1 9a ; fe 1 



TO,- H Trrrra U>4 H a,W 



so that 



r 9r r 00 fc 



^(r 2 ) = 2(r 2 + rr) = 8po + 0(r-^. (2.8) 
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Noting that T < can be chosen so that r(T ; xq, £°) is large and r(T ; xo,£,°) < 
0, we learn (|2.8|) implies (|2.6p . Taking a small neighborhood of expToiJ p (xo, 
the assertion on local uniformity follows. □ 

For any small open set W C T*M \ T_ and negatively large To < 0, we have 

exptff p (x , £°) = (r(t; x , 0(t; * 0) p(t; x , w(t; x , ?)) € T* Moo 

for (x ,£°) G W - and t < T . Then, identifying ^ M free , we can define for 
(x Q ,£°) € and t < T 

St(xo,€°) 

= cxp(—tH p 2) o exp ii? p (xq, £°) 

= (r(t; x , C°) ~ 2tp(t; x 0) 0(t; as , p(t; io, £°)> *o, G T*M free . 



Proposition 2.2 Let (x ,£°) G T*Af \ 71, f/ierc i/je firm* 
S_oo(zo,£°)= hm 5 t (i„,?°)erM frcc 

t — > — oo 

exists. Moreover, S-oo- T*M \ 71 — > T*Mf rcc is a local diffeomorphism. 
Proof. By the equations (|2.3|) . I|2.5p and the estimate (|2.7|) it follows that 



= —h lm u} m ——uj j Lu k + 0{r 1 ^h jk UjU} k ) 



2 Bh lm 

- - 1 W k ——uj l uj m uj k + Oir- 1 -*^ + y k uo 3 uo k )) 
r z oui 

= 0(r- 1 -»(l + h : > k Lj j u k )). 
Then by Lemma |2. II we obtain 

V k (0(t;x o , e Q )H-(i;xo,C K(i;xo,e°) < Cexp( / r(t; x Q , £ )- 1_/1 #) < C, 



which holds locally uniformly in the initial data. Thus uj(t;xQ,^°) is bounded 
as t — > — oo, and we have 

^-Mt; x , £°) - 2tp(t; xo, ?)) = 2a p + -a^ - 2tp = Oft*}- 1 -"), 
dt r 

so that there exists r_(xo,t; ) such that 

|r_ (x 0l (?) - (r(t; z , £°) - 2*p(t; x 0> < C<*>-' i . 
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Similarly for other variables there exist 6L(xo,£°), p (xo,^ ) and u> (xo,^°) 
such that 

\e^(x o ,e)-0(t;x o ,e)\ <C(t)-», 

|p-(xo,c )-p(t;^o,C )l<c(*)- 1 - /1 , 

\io'(x Q ,e)-^(t;x ,e)\ <c(t)-». 

All these estimates hold locally uniformly in (xo,£°), an d thus 
5-00= lim S t :T*M\T_^T*M hee 

t — > — oo 

is locally uniform convergence and S_oo is continuous. 
Next we prove the smoothness of S'-oo- Since 

S'-oo — Gxp(—tH p 'i) o S-oo ° exptHp for any tel (2-9) 

and exp£iJ p 2 and exptH p are diffcomorphisms on T*Mf roc and T*M, respec- 
tively, we have only to work on T*Moo C T*Mf lee , that is, it suffices to show 
that all the derivatives of St converges as t — ► — oo uniformly on any small open 
set W C T*M X with 

exp tH p {W) C T*Moo for t < 0. 

If we put 

(R(t),e(t),p(t),u(t)) =exp(-^ p2 )oexpti/ p ( a ;o,e ), (ar ,£°) € W, 
then it satisfies 

• <9X • cW . dK . dK . n . 

R = W 9 = ^ p = -tm> " = -&>> (2 ' 10) 

where 

K(t, R, 9, p, lo) = p 2 - p(R + 2tp, 9, p, 9) 

= - {R l Up)2 h^LjjLJk -a (R + 2tp, 9)p 2 

a{ {R + 2tp, 6) ai l k (R+ 2tp, 6) 

R + 2tp W {R + 2tpf 

Let us first show that all the derivatives of (R(t),9(t), p(t),w(t)) are bounded 
for t < 0. We have already shown the boundedness for (R(t),9(t), p(t),uj(t)). 
Now assume that {d^R{t), d p 6{t), d^p(t), d^u){t)), < I - 1 is bounded, where 
d 13 = d^ XQ for simplicity. Then, by differentiating (|2.10p . we learn for a G Z™ 
with |a| = Z 

= 0({t)- 1 ~' t {\(8 a R(t), d a 9{t), d a p(t),d a cu(t))\ + 1)), (2.11) 

at 



9 



and similarly for d a 6(t), d a p{t), d a u;(t). Thus, if we put 

m = \(d a R(t),d a e(t),d a p(t\d a u;(t))i 

then 

rO 

1-/*, 



ip(t)<C + CJ (sj-^ip^ds. 

The Gronwall's inequality applies and we get 

ifi(jb) < C. 

By induction all the derivatives of (R(t),0(t), p(t),u>(t)) are bounded, and thus, 
in turn, substituting this estimate to (|2.11[) and the counterparts for d a 8(t), 
d a p(t), d a cu(t), we can show that all the derivatives of St converges locally 
uniformly as t — ► — oo. 

Finally, in order to prove that S-oo is a local diffcomorphism, we need to 
estimate the first derivatives of St more precisely. Noting the boundedness of 
St(W), t < and the relation (|2.9|) , we observe it suffices to show that 5_oo 
is diffcomorphic on a small open set W = exipT H p (W) C T*M \ 71. For 
arbitrary large R\ > 0, by taking Tq < negatively large enough accordingly, 
we may assume 

S t (W) cT*((iJi,oo) x dM) \71. 

Thus similarly to (f2~TTj) we obtain for (i? , 6» , p , w°) € W 

d_ ( d(R(t),0(t),p(t)Mt)) 
dt\ d(R a ,9 0l p ,uj ) 



< 0((Rx+2tp) 



d(R(t),e(t), P (t)Mt)) 



and 



d(R ,9o,P°^°) 
d(R(0),6(p),p(0)M0)) 



d(R Q ,e ,p°,uj 

so that, if Ri is large, 

d(R(t),e(t), P (t)Mt)) 



d(R ,9o,p°,LO°) 



- 1 



= 1, 



< - for t < 0. 



It follows that the Jacobian of S-oo is invertible on W, and thus by 
is a local diffeomorphism. 



□ 
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3 Microlocal smoothing estimates 

Note that, if a boundary decomposition = (1. oo) x dM is fixed, the scaling 

T* Moo ->T* M x , (r,0,p,u) h-> (h~ x r,6,h~ x p,h- x u), ft e (0,1] 

is well-defined. Thus, if a € C^^M^), then a(hr,6,hp,huj), ft G (0,1] is 
well-defined. Choose ip G C°°(M) so that it does not depend on r near <9M, 
equals 1 on 

{(ft _1 r,0) € Moo; ft 6 (0,1], (r,0) G Tr(suppa)}, 

and has support sufficiently near the same set. Here ir: T*M — > M is the 
projection to the base space. If supp a is sufficiently small, we can choose ip and 
a coordinate neighborhood U C 9M such that supp ip C (1, oo) x U. Using the 
chart diffeomorphism k: (1,oo) x U ^> (l,oo) x (Jc M n , we have for v ETt 

K*[{K*g)~±{K^){K*a) w {hr,9,hD ri hDe){K.*^){K*g)^{K*v)] G 7Y, 

where k* and = (k*) -1 are the associated pull-back and the push-forward, 
respectively. In the following argument there always appears only one coordinate 
chart. Hence (l,oo) x U may be identified with (l,oo) x U C R", and k* and 
k* will not be written explicitly, if there is no confusion. 

Definition 3.1 Let (r , 9 , p° G r*Moo and u G S'(M). VKe write 

(r o ,0 o ,pV°)£ WF rl » 

if there are a G C^°(T*A/ 00 ) witfft a(r , O , p°, u°) ^ and ^ G C°°(M) as a&ove 
smc/i that 

\\g-i^a w (hr, 9, hD r , hD e )ipg*v\\n = 0{h°°) as ft | 0, 

(where we have omitted k). We call WF rh (n) the radially homogeneous wave 
front set of v. 

Remark Recall the semiclassical characterization of the (usual) wave front set, 
i.e., (xo,^°) WF(d) if and only if there are coordinate neighborhood V around 
x , a G C$°(T*V) with a(x ,S,°) ^ and ip € Cg°(V) with V = 1 on Tr(suppa) 
such that 

\\g-*iPa w (x, hD x )ip g *v\\ n = 0(ft°°) as ft J, 0. 

(See also the frequency set, Definition 13.41 ) The radially homogeneous wave 
front set differs from the wave front set only in the parameter in front of r in 
the polar coordinates. 

The "homogeneous" wave front set corresponds to a somewhat different set 
corresponding to the operators of the form 

a w (hr,9,hD r ,h 2 Dg). 

The homogeneous wave front set is used for the long-range case in [T3] and [B] . 
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Theorem 3.2 Let u e H, and t > 0. Suppose (x ,£°) € T*M \ 71, and 
denote its scattering data at t = — oo by (r_ , 6—, p~ , uj~). If 

(-2t Q p-,d^,p- 1 io-)^WF lh (u Q ), 

then we have 

(x 0i e)£WF(e- u ° H u ). 

Theorem l3.2l is a refinement of the results of [H] and [B] in the short-range case, 
and also they are refinements of the microfocal smoothing property of Craig- 
Kappeler-Strauss [jQ and Wunsch [21]. While [T] requires the initial data to 
decay in some cone in the configuration space in its assumption, the homoge- 
neous or radially homogeneous wave front set capture the decaying property in 
cones in the phase space, which microlocally weaken the assumption of [T]. 

Corollary 3.3 Let J be as in iTOj) . Then 

WF(e- ltH u ) \ 71 = WF(e- ltH J J*u ) \ 71 t > 0. 

Proof. Note that the operator J J* is the multiplication operator by j(x) 2 . Then 
we decompose 

e - ttH u = e- UH JJ*u + e- UH (l - JJ*)u . 
It is easy to see WF rh ((l - JJ*)u ) = 0, so that 

WF(e- M (l - JJ*)u ) \ 71 = t > 0. 

Then the assertion follows. □ 

The proof of Theorem 13.21 is analogous to [14j , and rather long. We first prove 
some lemmas needed later. 

Definition 3.4 Let Vh G TL be a vector dependent on h G (0,1] with \\vh\\ < 
1, and (xo,£°) £ T*M. We say that Vh is microlocally infinitely small near 
(xo,^°), if there exist a chart k: V — > V C E" around Xq, a G Cq°(T*V) wt/i 
a(a;o,^°) 7^ 0, and a cutoff function ip <E C^°(V^) wit/i ip = 1 on suppa suc/i i/iai 

|| 3 -iVa w (^^)V5^ft||w = 0(h°°). 
The frequency set FS(u/,,) is t/ie compliment in T*M of such (xo>£°) 's. 
Lemma 3.5 Let (x 0> (?) G T*M and veH. Then for any T Q G E 

(x ,e°)G WF(«) expT tf p (x ,£°) eFS(e-« ftr »^). (3.1) 

Proof. Let k: — ► V C E n be any chart around (xq,^ ). We assume To < 
and expiiJp(x ,£°) G V for t G [T ,0]. Suppose a G Cg°(T*V) has sufficientiy 
small support near (xo,£°), and put 

bo(x,Z;t) = aoexp(-tH p )(x,£) (x,£,t)eT*Vx [T ,0]. 
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Note that bo is the solution to the first transport equation: 



— 6 + {p, b } = 0, b (-,-,0) = a, 



where {•, •} is the Poisson bracket. Since suppa is small, we may assume 



Take tpo & Cq°(V) with ipo — 1 sufficiently near Ut e [r 0i o]7r(supp6(-, •;*)), and 



then, by the standard argument in the pseudodifferential calculus (see, e.g. [9]), 



-B (t, h) + ih[H, B {t, h)} = g(x)-*i/;i(x)rS(x, hD x -t, h)ipx{x)g{x)* , 



where r € S(h) is supported in supp&o modulo S(h°°), and tpi G Cq°(F) is 
any function with = 1 on supp^o- Here we have put S(h) — S(h, dx 2 + <i£ 2 ) 
in Hormander's notation in the chart we are concerned, and the dependence 
on the parameters t, h are supposed to be uniform. (This symbol class is not 
invariant under the coordinates change, but we do not have to take care of this 
fact, since we have fixed one chart.) Next we decompose ro = r' Q + r' ' so that 

suppr^, C supp& , r'JeSih 00 ), 

and consider the second transport equation: 

J^i + {p,&i} = 4 b 1 (-,-,0,h)=0. 
The solution is written as 



|J suppb (;-;t)cT*V. 



tepo,o] 



define 



B (t,h) = [g(x)-iMx)b^{x,hD x ;t)Mx)g(^] 2 , 




Let supp ipi be sufficiently near supp ipo and put 



Bi(t,/i) =g(x)-*^ 1 (x)bY(x,hD x ;t)Mx)g(x) 1 i 



then we obtain 



f t (B (t, h) + Bi (t, h)) + ih[H, B (t, h) + B 1 (t, h)} 
= g'^(x)tp 2 {x)ri(x, hD x ;t, h)ip 2 (x)g(x)* . 
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Here, similarly to the above, r\ S S(h 2 ) is supported in supp&o modulo S(h°°), 
and tp2 € C^{V) is any function with ip2 = 1 on suppi/^. Iterating this proce- 
dure and putting 

oo 

where ^ S Cq°(V) is any function with ip = 1 on U™ supp-0j, we obtain 

— ft) + ft)] = 0(h°°) in £(W), 

B(G,ft) = [ 5 (x)-3VoWa w (^,^x)^o(x)g(x)3] 2 . 

Also note that B(Tq, ft) is an ft-pseudodifferential operator with principal symbol 
aoexp(~T H p )(x^) 2 . Then 

(e- lhT ° H v, B(T , h)e~ lhT ° H v) n - (v, B(0, h)v) n 

= f T ° ^-{e- %htH v, B(t, h)e- MH v) H dt = 0{h°°), 
Jo dt 

and thus (|3.ip holds. 

For general To S R, divide the geodesic into a finite number of small segments 
each of which is contained in some coordinate neighborhood, and apply the 
above argument on each chart, then the lemma follows. □ 
Let uo € H, (xo, £°) E T*M\T- and to > be as in the assumption of Theorem 
By Lemma 13.51 it suffices to show 



ex V T H p (x ,e) i FS(e-^' +hT ^ H u ) 

for some To G K. Let us choose sufficiently large negative To *C so that we 
can work in a fixed coordinate neighborhood (1, oo) x U C M^, U C dM near 
the trajectory: 

ex.ptH p (x ,t°) £ (l,oo) x U for t < T . 

We put 

(r(t),6(t),p(t),Lj(t)) = exp(t + To)H p (xo,Si a ). 

Fix small 5' > and large C > 0, specified later. Take any So 6 (S'/2,S') and 
define 

<yj: (-oo,0] x T*M -> K 

by 

/ /, n / \r-r(t)\ \ ( \6-0(t)\ 
<p.(t,r,e,p,uj)=x[ ' u x' 



5<5o|* + T |y * ^-Clt + Tol-M^ 



5 - c|* + Td-A'-i ; A v^o - c\t + t \-k 

X1X2X3X4, 
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where 



x eC°°([0,+oo)), 



X(r) = 



1, ifr<l, 
0, ifr>2, 



(XT 



Lemma 3.6 Let small 5' and large C be given, and take T$ large enough accord- 
ingly. Then for any <5 € (S'/2, 5') is well-defined and satisfies the following: 



1. 



and 



<p(t, r, 6, p, ui) > o?i (-co, 0] x T*M, 
ip(t, r(t), 9(t), p(t), uj(t)) = 1 for all t < 0. 



2. The inequality 



D 

TTt 



ip(t, r, 6, p, w) < on (-oo,0] xT*M 



holds, where is the Lagrange derivative defined by 

D d d dp d dp d dp d dp d 

Di = m + {p, ' i = lTt + l^lfr + d^M~Tfrdp~~ ded^' 

3. The inequality 

\d k r d%d l p dtd^\ < c kalfjn (t)- n - k 

holds, that is, 

d?ip e S((t)- n , (t)- 2 dr 2 + d9 2 + dp 2 + duo 2 ) uniformly mt<0. 

Proof. 1. is obvious. 
2. Note that 



D 



D X2 . Dx3 . Dxa 

J^rX2X3Xi + Xl-^7-X3X4 + XlX2~^7-Xi + XlX2X3~^- 



Now let us compute the differentiations. We first get 



Dxi _ 1 
Dt 5S \t + T \ 



\r — r(t)\ r — r{t) (dp dp 



\t + T \ \r-r(t)\ \dp dp 



:(t) 



f \r-r(t)\ 
\h5o\t + To\ 



On supp(-^-X2X3X4) we have 



i< l :; r{t) [<2 

- 5S \t + T \ ~ 
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and 



dp _ dp 

Op d P [) 



^ASo + Odt + Tol- 1 ^). 



Thus, taking T larger if necessary, we obtain 



"X2X3X4 



< 



1 



Dt ~ 5S \t + T Q \ 

Similarly, by direct computations, 

Dx2 = 1 
Dt 5 -C\t + T \-» 

C^t + Tol- 1 -" 

D X 3 1 



(S + 0(\t + T \- 1 -n)x'X2X3X4 < 0. 



\0-6(t)\ 



+ 



6(t) (dp _ dp 



Dt 5 -C\t + T \-i-» 

C(1 + m)|* + T |- 2 -" 
1 



S - C\t + T \-» \0-6(t)\ \duj du 
\P-P(t)\ 



X 



P ~ P(t) ( dp dp 



DXi = 

Dt S -C\t + T \-» 

C l i\t + T \- 1 -' 1 
and, on supp(xiX2X3X4), 



5o-C\t + To\- x -n \p-p(t)\\dr dr 
\lj — u)(t)\ u) — u(t) (dp dp 



(*) 



5 -C\t + T \- 



\u-u(t)\ \d9 d9 



;(*) 



x , 



dp 


dp 


duo 


duj 


dp 


dp 


dr 


dr 


dp 


dp 


~bl ~ 


d6 



^Odt + Tol" 1 -"), 
<O(\t + T \- 2 -»), 
<O(|t + T |- 1 -' 1 ). 



Thus we have 

D X 2 



X1-^X3X4 < So _ C]t + To] -, [CM + Tor 1 -" + 0(1* + Tor 1 ^)] X1XX3X4, 

XiX2^X4 < So _ clt l Tol -i-» + ")l* + T or 2 ^ + 0(|t + T r 2 -")] X1X2XX4, 

[Cp\t + Tol" 1 "" + 0(|t + Tol" 1 "")] X1X2X3X'- 



^X4 . 
XlX2X3^^ < 



1 



Dt - S -C\t + T \-» 

Since C is large enough, we obtain the assertion. 
3. The estimates follow from direct computations. 



□ 
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We put 

b a {r,9,p,uj;t,h) = ip(hrH, r, 6, p, u) G 5(1, (h~H)- 2 dr 2 + d9 2 + dp 2 + dtu 2 ), 

and we restrict the parameter t G (— oo,0] to the interval [— io,0]. Let ipo G 
C°°(M 00 ) be such that it does not depend on r near dM, equals 1 on 

(J{tt(supp& (-, t,h)); te [-t ,0],he (0,1]} 
and has support sufhciently near the same set. We define 
B (t,h) = [g(r,eyiMr,e)b^(r,9,hD r ,hDg- 1 t,h)Mr,0)g(r,e)i] 2 on H. 

Lemma 3.7 There exists r G 5(1, (h~ 1 t)~ 2 dr 2 + d9 2 + dp 2 + duj 2 ) such that 

^-B {t,h)+i[H,B (t,h)] 
dt 

< g(r, 0)-*V>i(r, 6)r%{r, 9, hD r ,hD e ; t, fc)^i(r, %(r, 9)i, 

and that tq has support in supp&o modulo S(h°°). Here ipi G C 00 {M 00 ) is any 
function that does not depend on r near dM, equals 1 on suppf/'o? and has 
support sufficiently near the same set. 

Proof. We compute the principal symbol of 

j t B {t, h) + i[H, B (t, h)] = j t B (t, h) - i[A g , B (t, h)} + i[V, B (t, h)], 

and apply the sharp Garding inequality. Let ipi G C°°(M 00 ) be as in the 
assertion. By the standard argument in the pseudodifferential calculus we can 
write 

j t B (t,h)-i[A g ,B (t,h)} 

= g-hi[2bo(dtb + h-^p, b }) + r .i} w (r, 9, hD r , hD g ; t, h)^gi , 

where the remainder term ro,i G 5(1, (/i~-4)~ 2 dr 2 + d9 2 + dp 2 + dcu 2 ) has the 
support property 

suppro.i C suppbo mod 5(/i°°). 

It is also easy to see 

i[V, B (t, h)] = g(r, 9)-iMr, 0)r£ 2 (r, 9, hD r , hD g ; t, h)^{r, 9)g(r, 9)i , 
r .2 G S(h», (h- 1 ty 2 dr 2 + d9 2 + dp 2 + duj 2 ), 
suppr 0j 2 C supp6 mod 5(/i°°). 

Thus, by Lemma l3.61 the principal part is 2bo{dtbo + h^ 1 {p, bo}), and it satisfies 

2b (d t b + h^ip, b Q }) G S(h-\, (h-H)- 2 dr 2 + d9 2 + dp 2 + du 2 ), < 0. 
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By the sharp Carding inequality r ,3 G 5(1, (h 1 t) 2 dr 2 + d9 2 + dp 2 + du> 2 ) is 
found such that 

suppj'0,3 C supp&o mod S(h°°) 

and that 

2b (d t b + h~ x {p, b }) w (r, 9, hD r ,hD e ; t, h) < r™ 3 (r, 9, hD r , hD e ; t, h). 
Putting the remainder terms roj, j — 1, 2 together, we obtain 

^-B (t,h)+i[H,Bo(t,h)] 
dt 

< g(r, 9)~iMr, % >, 9, hD r , hD e ; t, ft)^(r, 9)g(r, 9)i , 

where r = r ,i + r ,2 + f*o,3- □ 
Let 5' , C and T be as in Lemma 13.61 Take an increasing sequence 

6' 

— < 5 < Si < 5 2 < ■ ■ ■ < 5 < 5', 

and set 

\r-r{hrH)\ \ ( \6-Q{hrH)\ 



bj(r,0,p,u;t,h) = x 



■ X 



bSjlh-H + To]) \5j - C\h-H + T \-v , 

\p- P {h-H)\ \ ( \u - uj{h-H)\ 



Jj-Clh-H + Tol- 1 -* J * \Sj -C\h-H + T \-i* 
Let tjjj G C°°(M 00 ) be such that it does not depend on r near dM and that 

^' = 1 on |J{7r(supp6 j (-,-, h)); t G [—to, 0], h G (0, 1]}, 

supp^j c y{7r(supp6j+i(-, M); ^ e [-*o,0],ftG (0,1]}. 

We set for j = 1, 2, . . . 

^.(t, ft) = 3 (r, 6>)-i^i(r, 9)C 3 h^ 3 -^ \t\bj(r, 9, hD r , hD B ;t, h^j^r, 9)g(r, 6)i. 

The constants Cj will be determined inductively. Since b\ is bounded from 
below by a positive constant on supp&o, there is large C\ > such that 

r (r,9,p,uj;t,h) < C\bx{r,9,p,u;t,h), 

where we have decomposed ro = r' + r' ' so that 

suppr-g C supp6 , r' ' G S(h°°). 

Now let us consider the operator -f^Bi + i[H, Bj], and iterate the argument in 
the proof of Lemma 13.71 Noting that 

j| G ^(r)- 1 ""^^) 8 , (r)- 2 dr 2 + d9 2 + (p;io}- 2 dp 2 + (p-^dcu 2 ), 

^ G S( (r)- 1 "^ (p; w) , (r) " 2 dr 2 + d9 2 + (p; uj) - 2 dp 2 + (p; to) ~ 2 duj 2 ) 
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and r = 0((h 1 t)) on supp&i, we can write as operators on L 2 (M. n ) through 
the chart diffeomorphism 



g^ef^-B, +i[H,B 1 ])g(r,6)-i 

7 

(r, 9, hD ri hDo; t, h), 



Ci|t|^-Ci6i+n,i 



where 



Since 



Ci\t\^ - C161 G S(l, (h-H)- 2 dr 2 + d9 2 + dp 2 + dcu 2 ), 

n,i e S^", (h-H)- 2 dr 2 + de 2 + dp 2 + du 2 ). 

suppri ; i C supp^! mod S(h°°). 



Ci|t|^-(C 1 6 1 -r' )<0 ) 



we can find by the sharp Garding inequality 

ri,2 G S(ft, (h-H)- 2 dr 2 + d6 2 + dp 2 + dcj 2 ), 
suppri i2 C suppbi mod S(h°°) 

such that 



(r, 9, hD r , hD g ; t, h) < rf 2 (r, 6, hD r , hD g ; t, h). 



Also here the inequality is in the operator sense on L 2 (M. n ). Combining with 
the remainder terms j-q and ri ; i, we obtain 



d_ 

dt 



Bi(t,/i) + i[JT, Bi(t,/i)] 



< </(r, 6)~*Mr, 0)rT(r, 9, hD r , hD g ; t, h)Mr, 0)g(r, 0)* 
- g(r, 0)-*Vo(r, 9)r%{r, 9, hD r , hD ;t, /i)Vo(r, %(r, 0)* 

as operators on H. Thus 

^(Bo(t, h) + ft)) + i[IT, B (t, ft) + 

< 5 (r, fl)-^i(r, 0K(r, 9, hD r , hDg- 1, /i)V>i(r, 6)g{r, 0)*. 



We repeat this procedure for j = 2,3,. 
the constants Cj > in the definition of Bj and 



By induction we finally obtain 



G S(hM-V t (h-H)- 2 dr 2 + d9 2 + dp 2 + duo 2 ), 
supprj C supp&j mod S(h°°) 
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such that 
d k 



- y £B j (t ) h) + i[H,Y t B j (t,h) 



3=0 3=0 

< g{r, 6)-fy k {r, 9)r^{r, 9, hD r , hD e ; t, h)^ k {r, 9)g{r, 0)i, k=l,2,.... 

Let ip G C°°{M 00 ) be such that it does not depend on r near dM, equals 1 
on U^g suppipj, and is supported in its sufficiently small neighborhood. 

Lemma 3.8 There exists an operator 

B{t, h) = g{r, d)-itp(r, 9)b w (r, 9, hD r , hD e \ t, h)^{r, 9)g{r, 
b G S(l, (h-H)- 2 dr 2 + d6 2 + dp 2 + duj 2 ) 

with the following properties: 
1. 

B(0,h) = B {0,h). 

2. For any e > we can choose b such that supp6(-, — {to + HTq), h) is 
contained in 

{{r,6,p,uj)eT*M 0o ; \r - 2/rV~*o| < ehr 1 ^ |0-0_|<e, 
\p - h~ l p-\ < eh- 1 , \u - h- x uj-\ < eh' 1 } 

modulo S{h°°). 

3. The Heisenberg derivative of B{t,h) satisfies 

5B{t, h) = j t B{t, h) + i[H, B(t, h)] < g{r, d)~h{r, 9)R{t, h)^{r, 9)g(r, 9)* 

on TL, where, through the chart diffeomorphism, R(t, h) is a bounded op- 
erator on L 2 (R") with sup_ to<t<0 || R{t, h)\\ = 0{h°°). 



Proof of Theorem \3.2l We have 

{e-^ ta+hT ^ H uv, B(0, h)e-^" +hT ^ H u )H 
= {u , B{-{t + hT ), h)u )n 

{e-^ +t ° +hT ^ H u , 5B(t, h)e-^ +t °+ hT ^ H u ) n dt 

(t +hTo) 

< 0{h°°). 

Thus the theorem is proved. □ 



20 



4 Proof of Main Theorem 

Letu e Hand t > 0. We suppose (x 0l £,°) € T*M\T-, and let (r_, 0_, p - , w _ ) 
be the scattering data at < = — oo. As in the previous section, let To *C be 
large enough that we can choose a coordinate neighborhood (1, oo) X U C Moo, 
[/ C dM such that 

7to expt/f p (a: ,f ) 6 (2,oo) x U, < < T , 

and, moreover, 

7to S t > oexpiiJ p (a;o,C ) € (2,oo) X U, t< T , i' < 0. 

By Lemma 13.51 and its analogue on Aff rco it suffices to show that 

e^pT H p (xo,e) G FS(e- l ( to+,lTo ) H W o) 

(r_ + 2T o p-,0_,p-,ar) G FS( e - l(to+/lTo)ffo J*u ). 

Let a € Co°(T*((2, oo) x {/)) be a symbol supported near expT i? p (a;o, £°), and 
y € Cq°((2, oo)xU) be a function such that it is equal to 1 on Ut'<o^°St' (supp a) 
and is supported in its small neighborhood. We define the operator A(h) by 

A(h) = K*(K*g)~4 (K*ip)(K*a) vr (r,0,hD r ,hDe)(K*ip)(K,*g)i K* on H, 

where k: (1, oo) X U — ► (1, oo) X U C K" is a chart diffeomorphism, and will be 
omitted, as in the previous section. Then, putting t' Q = t + hT for simplicity, 
we have 

\\A{h)e~< H u \\ 2 H = {e-< H °J*u ,B{t' ,h)e-< H °J*u ) bee 

+ (e-^ H u , A(h) 2 e-< H (1 - JJ*)u )h (4-1) 
+ {A(h) 2 e-< H (1 - JJ*)u ,e-< H JJ*u )n, 

where we set 

B(t\ h) = e- lt ' Ho J*e lt ' H A{hfe- lt ' H Je lt ' Ha on H free 

for t' > 0. We first show that the last two terms in the right-hand side of (|4.1|1 
can be ignored: 

Lemma 4.1 For any T\ > 

\\A(h) 2 e- u ' H (l - JJ*)\\ H ^n = 0(h°°) uniformly in t' e [0,Ti]. 

Proof. Let A(t',h) denote the operator B(—t',h) constructed in the previous 
section with < t' < T x . Then 

(e~ lt ' H (1 - JJ*)v, 1(0, h)e- lt ' H (1 - JJ» H 

< ((1- JJ>,i(t',ft)(l- JJ*)v) H + T 1 sup ||i?(-t',ft)||||(l- JJ>|&. 

0<*'<Ti 
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Noting the supporting properties of operators A(t' , h) and 1 — J J*, we conclude 
that the right-hand side is 0{h°°) uniformly in t' G [0,T X ]. □ 

In particular, the second and the third terms in the right-hand side of (|4.1[) are 
0(h°°). Thus it suffices to show that B(t' , h) is an /i-pseudodifferential operator 
with principal symbol (ao S~].i //J 2 : since then the principal symbol approaches 

(a o Silo) 2 as h > tends to 0. 

Lemma 4.2 Let H be the self-adjoint differential operator on 7if rcc defined by 

H = J*HJ + (1 - J*J)H Q {1 - J* J). 
Then for any T\ > we have 

B(t', h) = e -it'H°jt'H( J* A ^2 j} e -it'Hjt'Ho + QQfo) m £(ft ft . oc ) 

uniformly in t' € [0, Ti]. 

Proof. Similarly to Lemma l4.ll we learn by duality that 

||(1- JJ*)e it ' H A(h)\\ H ^ H = 0(h°°) 
uniformly in t' £ [0,T±]. Since 

( J*H ~ HJ*)e u ' H A{h) = (J*H + (1 - J*J)H J*)(l - JJ*)e lt ' H A(h), 
we obtain 

\\e~ lt ' A (J*H - HJ*)e u ' H A(h)\\ H ^n = 0{h°°). 
Hence, by integrating it with respect to if, we have 

|| J*e lt ' H A{h) - e lt J* A(h)\\n^ H = 0(h°°). 

a 

Thus we have reduced the proof to the study of the behavior of 
B(t',h) = e- i *' ffo e lt '^(JM(^) 2 J)e" l *'^e i *' ffo . 

We note 

A_ e -it'H ^it'H = _i( e -i?H ( Ho _ jjyt'H j e -it'H 0e it'H 

= -iL{t')e- lt ' Ho e lt '" . 
If we write H - H = w w (r, 9, D r , D e ), then 

L(t') = w w (r-2t'D r ,6,D r ,Dg) 
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with no error terms by virtue of the Weyl calculus. Now B(t' , h) satisfies the 
Heisenberg equation: 

j t B(t',h)+i{L(t'),B(t',h)]=0, 
B(0,h) = [h(e)-i l fi{r,e)a w (r,6,hD r ,hDe)(p(r,6)h(9) 1 ^} 2 . 
Let us construct an asymptotic solution to this equation. We put 

b (r,6,p,uj;t',h) = (a o Sz\, /h {r, 9, p, uj)) 2 
for t' e [0,Ti], h e (0,1]. Note 

b (;;;-,t\h)eC^((2 7 ^) X U) , b G S(l), 

<p=l on |J{7r(supp6o(-,-,-,-;i',/i)); t' e [0,7\], fee (0,1]}. 

Since we are considering operators on Mf ree , we may use chart neighborhoods 
of the form R x U, U C in what follows. Choose a cutoff function ip n e 
C°°(R x [/) which does not depend on r, equals 1 on supp<p and is supported 
in its sufficiently small neighborhood, and define 

B (t', h) = h(0)-*MO)l>o(r, 9, hD ri hD 9 ;t' , h)M^M^ k - 

Lemma 4.3 There is r € ^((/i -1 ^} -1- ^) which is supported in supp b modulo 
S(h°°) such that 

^B (t',h) + i[L(t'),B (t',h)] 

= /i(0)"3 Vi(0)r^(r, 0, fc£> r) hDg; t', h)^(6)h{0) * . 

.ffere G C°°(]R x U) is any function that does not depend on r and equals 1 
on supp tp . 

Proof. Differentiating 

bo(S- t >/h(r, d, P, w); t' , h) = a(r, 9, p, uj) 2 
with respect to t', we obtain 
d 

—b + h{k, b } = 0, b (r, 9, p, uj; 0, h) = a(r, 9, p, uj) 2 , 

where k is the principal part of w(r — 2h~ 1 t'p, 9, h~ 1 p, /i _1 o>): 

k{r, 9, p, uj; t' , h) = j(r - 2h' 1 t' p) 2 [h- 2 p 2 - p(r - 2h~H' p, 9, h' 1 p, h' 1 ^)}. 
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Since r, p,u> — 0(1) and p < c < for (r, p, u>, t' , /i) € supp 60, we have 

\d J r d^d k p d^w(r^2h- 1 t'p,9,h- 1 p,h- 1 uj)\ < ChT 2 (rV)- 1 "" on su PP 6 , 

(4.2) 

and thus, if we apply the asymptotic expansion formula, 
-^B (t',h)+i[L(t'),B (t',h)} 

N-l 

, r -1 W - 

= h{9)-^M0)[Y. h j r 0d +R N \ {r,6,hD r ,hD 8 ;t',h)ipi(9)h(9)* 
3=0 

with 

Therefore it suffices to confirm that the remainder term Rn has an appropriate 
decaying property in h, since, then, the required ro will be obtained by taking 
the asymptotic sum of roj. However, since the estimate (|4.2j) does hold on the 
whole R 2 ™, we can not estimate it directly. Let us rewrite 



i[L(t'),B (t',h)] = ie- ltHo [H - H,e itH °B (t',h)e- ltHo }e ltHo 



°B (t', h)e~ ltH ° = h(6)-iM6)b™(r + 2t'D r , 9, hD r , hD g - 1', h)i/> (9)h(0)* , 



and compute the remainder term of h(9)i[H - H, e itH °B Q (t', h)e~ itH °]h(9yi , 
instead. The 7V-th remainder term is given by 

n N [w(r, 9, h~ x p, h~ 1 uj)b {) (r l + 2h' 1 t'p 1 , 6 1: Pl , un; t', h) 
- b (r + 2h-H'p, 9 : p : w; t', h)w(n, B x , h^pi, /i~Vl)] 

(ri,e»i,pi,wi) = (r,0,p,w) 

(4.3) 

modulo S^/i 00 ) (see, e.g., [5]). Here 

N ~ Jo (N-l)\ 

r ih ~\ N 
■ |y(dr<9 Pl + dgd^ - d p d ri - d^dejj , 

which is defined as a Fourier multiplier by 

J! ^^e^^+ 6a ^- a ^dr{^(m + §h ph - ^" 1 )} JV . 

Write (|4.3p in the integral form with respect to (r, 9, p, u>, T\, 9i, pi, u>i) and their 
conjugate variables (r, 6, p, cj, fi, 61, pi, Cj\), and integrate it by parts to make it 
intcgrable. By the supporting property of bo we obtain, for example, 

d J r d%d k p d%{r + 2h-H' p,0, p,u;t' ,h) G S((r) k ,dr 2 + d9 2 + h- 2 t' 2 dp 2 + dcu 2 ), 
w(r, 9, h~ x p, h^uj) G S(h~ 2 (p;uj) 2 , (r)- 2 dr 2 +d9 2 + (p)- 2 dp 2 + (uj)- 2 du 2 ), 
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and we learn that (|4.3[) belongs to 

S(h N ~ M , dr 2 + d9 2 + h- 2 dp 2 + dcj 2 ), 

where M is chosen independent of N. Then it is easy to see that the asymptotic 
expansion is verified. □ 

Now we solve the second transport equation: 
d 

— hi + h{k,bi} = -ro, h(r, 9, p, u; 0, h) = 0. 
The solution is given by 

h(r, 9, p,w;t',h) = / r (S- s / h o Szl,/ h (r,0,p,io);s,h)ds, 
Jo 

and it satisfies by Lemma [4731 

bi£S(h), supp b\ C supp b Q mod S(h°°). 

Similarly to Lemma [4~3l there is r-y £ SQifi^t') -1 ^^) such that it has support 
in supp&o modulo S(h°°) and that 

^(Bo(f, h) + Bitf, h)) + t[L(t r ), B (t>, h) + B^t', h)} 
= h(9)-iiP2{6)rJ(r, 9, hD r , hD ;t', h)^ 2 {9)h{9)* , 

where 

Bx(t', h) = h(6)-iipi(6)bY(r, 9, hD r , hD e ;f, h)ipi(9)h(6)i 

and tp2 £ C°°(M x U) is any function that does not depend on r and equals 1 
on supp ipi . 

Repeating this procedure, we can inductively construct the symbols bj and 
rj. Let ip £ C°°(R x U) be independent of r and equal to 1 on U!*i supp , 
and set 

oo 

B(t', h) = h(9)-^ij(9)b w (r, 9, hD r , hD e ; t', h)ip{9)h{9)? , b-J^bj, 

3=0 

then we have 

^B(t',h)+i[L(t'),B(t',h)]=0(h™), 
B(0,h) = [h{9)-iip(r,9)a w (r,9,hD r ,hD e )ip(r,9)h(9) 1 i} 2 . 
This implies 

B(t', h) = B(t', h) + 0{h°°) = B(t', h) + 0{h°°). 
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Substituting this to (|4.ip , we obtain 

\\A(h)e-< H u \\ 2 n = (e-^ruo^it^h^-^ruo^ + Oih™). 

B(t' , h) is an ft-pseudodifferential operator with principal symbol (a o SZ t i /ft) 2 - 
Assume (r_ , 6>_, p~, u~) £ WF(e~ zt ° H ° J*u ), or equivalently 

(r_ +2T o p-,0_,p-,cj-) £ FS(e^ to+,lTo ^ J* W o). 

If the symbol a is supported sufficiently near expToH p (xo, £°), then the full 

symbol of B(to + hTo,h) is supported near (r_ + 2Top~ , 8-, p~ , lo~) for small 
ft, > 0, so that 

(g-iCto+hToJffo j* MQ) ^( tQ + hTo: h)e-^ to+hTo]Ha J*u ) hcc = 0{h°°). 
Thus we obtain 

\\A{h)e^+ hT ^ H u \\ n = 0{h°°), 

which implies expT H p (x , £°) ^ FS(e~ ,;(to+/lTo)ff u ). Hence it follows that 
(x ,£°) (£WF{e- UoH u ). 

Conversely, assume (xo,£°) ^ WF(e~ ltoH uo). We choose a symbol a > 
supported in a small neighborhood of expToH p (xo, £°) and construct the 

operator B(to + ft To, ft) as above accordingly. Then by the assumption we have 

{e- l{t0+hTa)Ho ru ,B{tv + hT Q ,h)e' l(to+hTo)Ho J*u ) ircc = 0(h°°). 

The principal symbol of B(to + ftTo, ft) is positive near (r_ + 2Top~, 8-, p~ , uj~) 
for small ft > 0, and thus it follows that 

(r- + 2T Q p-,9_,p-,ij-) i FS(e-« to+hT ^ Ho J*u ), 

or equivalently (r_, u~) £ WF(er ttaHa J*u ). The proof of Theorem O 

is complete. 
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